Spectral densities of Kondo impurities in nanoscopic systems 
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We present results for the spectral properties of Kondo impurities in nanoscopic systems. Using 
Wilson's renormalization group we analyze the frequency and temperature dependence of the im- 
purity spectral density pi(u,T) for impurities in small systems that are either isolated or in contact 
with a reservoir. We have performed a detailed analysis of the structure of pi(ui,T) for u) around 
the Fermi energy for different Fermi energies relative to the intrinsic structure of the local density of 
states. We show how the electron confinement energy scales introduce new features in the frequency 
and temperature dependence of the impurity spectral properties. 

PACS numbers: 72.15.Qm, 73.22.-f 
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I. INTRODUCTION 

The physics at the nanoscale has emerged as one 
of the most active areas in condensed matter physics. 
This new field includes the study of small metallic and 
superconducting islands, quantum dots and nanostruc- 
turated semiconductors, nanotubes and quantum corrals, 
nanoelactronics and nano-characterization among other 
things B 

The advances in nanotcchnologies revived the interest 
in the Kondo effect ,□ one of the paradigms of strongly 
correlated systems. On one hand, Scanning Tunneling 
Microscopy (STM) allowed the direct measurement of lo- 
cal spectroscopic gipperties of Kondo impurities on no- 
ble metal surfacesau and in nanoscopic systems. On the 
other hand, it has been shown that single electron transis- 
tors and single walled carbon nanotubes weakly coupled 
to contacts may behave as Kondo impurities generating 
new alternatives to study the phenomenal] 

In its simplest form, the Kondo problem is a single 
quantum spin interacting with an ideal electron gas. An 
antiferromagnetic coupling between the impurity and the 
free electron spins gives rise to an anomalous scattering 
at the Fermi energy leading to a large impurity contri- 
bution to the resistivity. Simultaneously the impurity 
spin is screened by the conduction electrons and the mag- 
netic susceptibility saturates at low temperature. There 
is a characteristic temperature Tk that separates the 
low temperature from the high temperature regimes. At 
T Tk, the impurity spin is essentially free and the 
problem can be treated by perturbations in a dimension- 
less coupling constant A. At T <§; Tk the impurity spin is 
screened forming a singlet complex with the conduction 
electrons and the system is described by an infinite effec- 
tive coupling. The crossover regime with T ~ Tk is more 
difficult to describe and the best treatment corresponds 
to the numerical renormalization group (NRG) as done 
by Wilson.Q 

The characteristic Kondo temperature is given by Tk 
~ D\f\e~ 1 l x where 2D is the free electron bandwidth. 
Associated to this energy scale there is a characteris- 
tic length scale known as the Kondo screening length 



£x = hvp/ UbTk where vf is the Fermi velocity. The 
physical meaning of the screening length is that, in 
the low temperature regime where the impurity spin is 
screened, the antiferromagnetic correlations between the 
impurity and conduction electron spins extend up to a 
distance of the order of £ k ■ 

The problem of Kondo impurities in nanoscopic sys- 
tems for different experimental realizations has been 
the subject of many recent theoretical and experimental 
works.Q When a Kondo system, either an atomic impurity 
or an artificial atom like a quantum dot (QD), is embed- 
ded in a small system of volume L d where d is the spacial 
dimension, the length-scale L should be compared with 
the characteristic Kondo length For L < finite- 
size effects arc expected to be important. The condition 
L ~ £,k is equivalent to UbTk ~ A where the energy 
A gives the average level spacing of the finite system. 
For a finite system the characteristic energy A acts as a 
low energy cutoff for the charge and spin excitations and 
consequently it modifies the low temperature behavior of 
the system. 

The ground state properties of a Kondo impurity in 
a small system have been addressed bM-,«unumber of au- 
thors using different approximations.l§cHl3 The thermo- 
dynamic properties and the effect of coupling the system 
to a macroscopic reservoir have recently-been studied us- 
ing the Wilson's renormalization group.B In this work we 
extend renormalization group calculations to evaluate the 
impurity spectral density and analyze how the new en- 
ergy or length scales introduced by the finite size affects 
the Kondo resonance at the Fermi energy. We also study 
the temperature dependence of the low energy spectral 
density. 

The rest of the paper is organized as follows: in section 
II we present the model and describe how the numerical 
renormalization group is adapted to our case. We then 
recap the most relevant thermodynamic properties. Sec- 
tion III contains the spectral properties of the impurity 
for frequencies close to the Fermi energy. After present- 
ing the general formulation we show results for finite sys- 
tems, systems in contact with a reservoir and the tem- 
perature dependence of the Kondo resonances. Finally 
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section IV includes a summary and discussion. 

II. MODEL AND THERMODYNAMIC 
PROPERTIES 

In this section we present the model for an impurity in 
a nanoscopic system coupled with a macroscopic reservoir 
and briefly discuss how Wilson's renormalization group 
is adapted to this problem. Then we review the thermo- 
dynamic properties of the system. 

A. Model Hamiltonian and Wilson's 
Renormalization Group 

Our starting point is the Anderson model for magnetic 
impurities pwith a Hamiltonian which in the usual nota- 
tion reads 

Ham = ^ e d d\d a + U d\d 1 d\d l + e v c\ a c urj 

+ J2( V » c ^ d ° + F *4<W) - ViBS iz , (1) 

where the operator <ij. creates an electron with spin a at 
the impurity orbital with energy e ( i and Coulomb repul- 
sion U, and c\, a creates an electron in an extended state 
with quantum numbers v and a and energy e v . The last 
term represents the effect of an external magnetic field 
along the z-direction coupled to the impurity spin S,;. 
Hereafter we will use D = 1 as our unit of energy. 

In this notation, the nanostructure of the system is hid- 
den in the structure of the one-electron extended states 
with energies e„ and wavefunctions ip u (r). In equation 
(|I|) the hybridization matrix elements are taken propor- 
tional to the extended state wavefunctions at the impu- 
rity position, i.e. V„ = VoVv(O) where the impurity posi- 
tion is defined as the origin of coordinates. We consider a 
simple structure consisting of a spherical metallic cluster 
of radius R c with the impurity at the center. The cluster 
is embedded in a bulk material with which it is weakly 
coupled through a large surface barrier. The Hamiltonian 
can be put in the form Ham = H AM + W(r,R c ) where 
the first term is the Anderson model Hamiltonian for an 
impurity in an infinite homogeneous host and the last 
term is a spherically symmetric potential barrier placed 
at a distance R c from the impurity. Figure |l|(a) illus- 
trates the configuration described by the model. For an 
infinite impenetrable barrier, the central cluster is de- 
coupled from the macroscopic reservoir and the model 
describes an impurity in a small sample. In this situation, 
the extended states that are coupled to the impurity, are 
confined in the central cluster and their energy spectrum 
is a discrete spectrum with a mean energy level sepa- 
ration given by the characteristic energy A. For a finite 
barrier these states are hybridized with the continuum ac- 
quiring a finite lifetime, the local density of states inside 




FIG. 1: A sketch of the central grain embedded in a metal 
(a). Inside the grain the characteristic energy level spacing 
A is shown together with the impurity levels at Ed and Ed + 
U. In (b) the linear chain obtained after Wilson's canonical 
transformation. 



the central cluster then presents resonances separated by 
the energy A and with widths 7 that are determined by 
the barrier. The model then incorporates the new energy 
scales A and 7, that may drastically change the impurity 
behavior. 

This geometrical structure is particularly appropriate 
to use the numerical renormalization group approach de- 
veloped by Wilson. Wilson introduced a logarithmic dis- 
cretization of the energy of the conduction electrons, di- 
viding the band in a series of energy intervals with ex- 
ponentially decreasing width. By means of a canonical 
transformation the resulting Hamiltonian can be mapped 
into a linear chain with variable hoppings. Each site rep- 
resenting an orbital surrounding the impurity with an 
associated energy (or length) scale. Wilson proposed 
to solve the linear chain by iterative perturbation. A 
truncated chain with N sites, described by an effective 
Hamiltonian Hn, gives the correct physics on an energy 
scale (Dat. A renormalization group transformation cor- 
responds to adding a site to the chain and relates the 
Hamiltonians describing successive lower energy scales. 
This leads to a systematic way of calculating the ther- 
modynamic properties at successive lower temperatures 
and the spectroscopic properties at successive lower fre- 
quencies. 

In our case, the Hamiltonian H\ M is rewritten in Wil- 
son's basis as schematically shown in Fig.^J(b). The po- 
tential barrier is described including a higher diagonal 
energy to the orbital centered around R c . Due to the 
structure of Wilson's basis wavefunctions, a given po- 
tential barrier leads to a diagonal energy that decreases 
as R c increases. Alternatively the barrier can be simu- 
lated with a smaller hopping matrix element in the region 
where W(r, R c ) is different from zero. We adopted this 
last description reducing one hopping term by a factor a. 
In what follows it is assumed that the band of extended 
states is half filled - the Fermi energy is set equal to zero 
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the broadening of the central peak for odd N c or the tails 
of the states at ±A/2 for even N c . The local densities 
of states as obtained with the NRG are also shown in 
Fig. [| For a^O the number of electrons in the central 
cluster is no longer a good quantum number and in what 
follow we refer to the situations of Figs. ||(a) and ||(b) as 
the Fermi energy being at a resonance (at-resonance) or 
between two resonances (off-resonance) respectively. 

As we show below, this structure of the local density 
of states determines the thermodynamic and the spectral 
properties of the impurity. 



FIG. 2: Unperturbed local density of extended states at the 
impurity coordinate (site 1) for a cluster with a — (vertical 
bars) and a — 0.05 (lines). In the later case the spectral 
density has been evaluated using a small imaginary part in 
the frequency (a) iV c = 19 and (b) N c = 20. 



- and td = —U/2. This guarantees that the electron-hole 
symmetry of the problem is preserved. 

The properties of a Kondo impurity depend on the lo- 
cal density of states at the impurity coordinate that in 
Wilson's representation is the local density of states at 
site 1. We end this section with a brief discussion on the 
effect of the confining potential W(r, R c ) on the local den- 
sity of states close to the Fermi energy. For a given N the 
spectrum of Hn depends on the parity of N. In the ab- 
sence of impurity, a one particle state at zero energy ex- 
ists only for odd N as schematically shown in Fig. g. The 
spectrum, as a function of N, alternates between those 
of Figs. g(a) and ||(b). The mean energy separation, 
between the one-electron states, is given by the charac- 
teristic energy A that decreases as N increases. This is 
the spectrum of a finite system (JV = N c ) described by 
an infinite barrier. The other one electron states, not 
shown in Fig. ||, arc not at energies nA or (n + 1/2) A 
with integer n. This is due to the logarithmic discretiza- 
tion of the band. How the NRG can be-used to analyze 
finite size effects was shown by NozieresEI. To show that 
the logarithmic discretization correctly describes the low 
energy spectrum of a impurity in finite systems, in the 
next section, we compare the exacts results on a finite 
system with those obtained with the NRG. The alternat- 
ing spectrum with a fixed Fermi energy corresponds to 
an alternating parity in the number of electrons, in fact 
for a fixed electron density, as the radius R c of the clus- 
ter increases, the number of electrons alternates between 
even and odd. We recall that due to the symmetry of 
the system only the sector of the Hilbert space with s- 
wave electrons is considered here. With a finite barrier 
(a ^ 0) the discrete states of the cluster are hybridized 
with the continuum of the host metal and become reso- 
nances. The local density of states then becomes a con- 
tinuum but retains some structure characterized by the 
energy A. In the NRG approach, Hamiltonians Hm with 
N > N c accumulate states at low energies representing 



B. Thermodynamic Properties 

Here we briefly recap the thermodynamic properties of 
the model. As stated above, for a = the at-resonancc 
situation corresponds to an odd number of extended elec- 
trons. An impurity in the Kondo limit contributes with 
an extra electron and the ground state is a singlet indi- 
cated as |0). The expectation value of the impurity spin 
(Q|(Sjz|Q) is zero reflecting the complete screening of the 
impurity spin and the zero temperature susceptibility is 
finite. For the off-resonance situation the ground state of 
the isolated cluster with a Kondo impurity is a Kramers 
spin-1/2 doublet, | -ff) and | JJ-). The expectation value 
(fl \Siz\ 1T) = —(4 \Siz\ -U-) i s different from zero and, 
in the low temperature limit the impurity susceptibility 
diverges as x = Mi (ft I Si* I 1\} 2 /k B T. 

For a finite barrier the susceptibility always satu- 
rates, however the at-resonance and off-resonance situa- 
tions give rise to very different temperature dependences. 
We-calculated the impurity magnetic susceptibility given 
byN 



k B T X = I 

V 

2k B Tj2\{»\Si z \0\ 



2 *v 



E% — E v 



(2) 



where the summation is done over the low energy states 
\v) with energies E v and P„ = ex${—E v /k B T)/Z. The 
matrix elements (HSi z |£) have to be evaluated in a recur- 
sive way at each renormalization step. The susceptibility 
reflects the thermodynamic properties of the system and 
we use the effective magnetic moment /i 2 = k B T\ as an 
indicator of the degree of screening of the impurity spin. 

In absence of barrier - corresponding to the infinite 
homogeneous system - the characteristic energy scale is 
the Kondo temperature indicated as k B T^ . A finite 
barrier at R c - enclosing N c shells- introduces the new 
energy scale A = DA~ Nc ^ 2 , where A = 2 is Wilson's 
discretization parameter. For k B T£? 3> A the fine struc- 
ture (on the scale of T£?) of the density of states does 
not change the properties of the system. Conversely, for 
k^T^ 3 ~ A new confinement induced regimes are ob- 
served: for the system at-resonance, as the temperature 
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k B T/T K 

FIG. 3: The magnetic susceptibility ksTx versus temper- 
ature. The thick line corresponds to an impurity in a bulk 
material (a = 1), lines and symbols to a barrier enclosing 
N c — 19 and 20 shells (open and full symbols respectively). 
Triangles (circles) correspond to a — 0.35 (0.6). The other 
parameters are: e d = —0.5, U = 1.0, and V = 0.2. 




FIG. 4: The effective Kondo temperature, that describes the 
low temperature behavior of the susceptibility, versus the bar- 
rier parameter a for the at-resonance (upper curves) and the 
off-resonance (lower curves) cases. Open symbols were ob- 
tained using Eq. (^), full symbols using the scaling shown in 
the inset. The parameters are the same as in Fig. H. 



is lowered, there is a rapid decrease in the magnetic mo- 
ment p? when ksT ~ A; for the system off-resonance, as 
ksT approaches A the magnetic moment saturates lead- 
ing to a plateau in the temperature dependence of p, 2 , 
only at lower temperatures the screening is completed 
(see Fig. ||). This plateau can be interpreted as the 
behavior of an isolated cluster. Only at very low temper- 
atures the coupling with the host states becomes relevant 
and the complete screening can occur. The condition for 
the existence of these new regimes fc^T^? ~ A can be 
put in the form ~ R c where f|? = Hvf /IceT^? is 
the Kondo screening length of the infinite system. In 
other words, only if the Kondo screening length is of the 
order or larger than the system size, the confinement ef- 
fects become evident in the thermodynamic properties. 
Although the behavior of the system is not universal, at 
very low temperatures the susceptibility can be used to 
define an effective Kondo temperature. In fact the low 
temperature tail of ksTx can be scaled to define the 
effective energy scale Tjj/' for the low energy spin exci- 
tations. We have done this using two different criteria, 
one is simply to use Wilson's result: 

X (T = 0)= 0.103-^, (3) 

to relate the low temperature susceptibility to an ef- 
fective Tgff. The other alternative is to plot /i 2 vs. 
T/Tk* fitting Tj[ to have a good scaling at low tem- 
peratures. The effective Kondo temperatures describing 
the low temperature behavior are shown in Fig. |J as 
a function of the barrier height for the at-resonance and 
off-resonance situations. For a > 0.2 the two criteria give 
the same results, for smaller values of a the universal be- 
havior is obtained only at extremely low temperatures 
where numerical errors become important. 



III. SPECTRAL DENSITIES OF IMPURITIES 
IN NANOSCOPIC SYSTEMS 

Using the standard definitions and notation, the impu- 
rity Green's function can be written, using the Lehmann 
representation, as: 

1 _ -E x /k B T , -Ey/kuT 

aw -g^ElPIWI' u _ ( l_ Ex) ■ 

(4) 

where Z(T) is the grand partition function 

Z{T) = Y J ^ Ex/kBT - (5) 

A 

The corresponding impurity spectral density is 
Pi(J {uj,T) = --[G ff (^ + ^0 + ,T)-G^-^0 + ,T)] 

= zlf) Ei M ^i 2 (^ /fcflT + e - £WfcBT ) x 

S(u - (E y - Ex)), (6) 

with Mx,x> — (Aider | A'). In the absence of an external 
magnetic field, the impurity spectral density is spin in- 
dependent and from here on we drop the spin index in 
Pia(w,T). 

We now briefly discuss the application of the NRG to 
evaluate these quantities. The zero temperature limit of 
expression (^) can be put in the form: 

Pi(u,T = Q) = -L^ |M A , | 2 <5( W + £ A ) 

^ ' A,0 

+ \M ^\ 2 S(to - E x ), (7) 

where the subindex indicates the ground state, the sum- 
mation is over all excitations A with energies E\ and on 
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the components of the ground state in the case of a de- 
generate ground state, Z(0) is the zero temperature par- 
tition function that gives the ground state degeneracy. 
We recall that in the NRG the chemical potential is set 
at zero and all many body energies E\ are measured from 
the ground state energy. 

In order to evaluate the impurity spectral density of 
the real system using the results of the NRG, consider the 
spectral density pf{u,T) corresponding to Hamiltonian 
i/jv- As mentioned above, the RG strategy of an iterative 
diagonalization of a sequence of Hamiltonians H]\f with 
N = 0, 1, ... is based on the fact that, on an energy scale 
u>n, the spectrum of Hn is representative of the spectrum 
of the infinite Hamiltonian £L Hence, for uj rs ujn it is a 
good approximation to takdi-J 

Pi (uj,T = 0)=pf(w,T = 0). 

A typical choice is to consider for the N th iteration the 
energy range A _JV / 2 < w < 2Ar N f 2 , that is to consider 
all the states of Hm with energies E^ in this range and 
all the corresponding matrix elements M^ to compute 
the spectral density as 



Pi(u,T = 0) 



1 



mo) 

rJV |2 



+\M^(u>-E^) 



(8) 



Since for the infinite systems the energy spectrum con- 
sists of a continuum, at each energy scale lon, the discrete 
spectra are usually smoothed by replacing the delta func- 
tion S(u) — E^) by a smooth distribution P N (ui — E^). At 
each N we use a Gaussian logarithmic distribution with a 
width that decreases as the characteristic energy scale of 
Hn- It is then clear that at high energies, say the atomic 
energy of the impurity level Ed, the NRG cannot resolve 
any detailed structure, only at low frequencies there is 
enough resolution to see, for example, confinement ef- 
fects. We analyze in detail the impurity spectral density 
around the Fermi energy which in fact contains all the 
physics at an energy scale relevant for the Kondo effect 
and for the electron confinement in a nanoscopic sample. 

Before presenting the renormalization group results we 
show, as a reference calculation, some exact results in 
finite systems. 



A. Exact results in finite systems 

Here we present the zero temperature results for an 
impurity in a finite system consisting of a linear chain of 
N equivalent sites with the impurity at one end. In this 
case the energy spectrum of extended states e v and the 
hybridization matrix elements V v of Hamiltonian (1) are 
given by: 



and 



-2tcos(W/(JV+ 1)), 



V v = V simW/(iV + 1)), 




1 o 1 
co/A 



FIG. 5: Impurity spectral density pi(uj,T = 0). Upper 
panels are exact results for a linear chain with N = 8 
(Vo = O.Ut, e d = E F - 0.5t, U = l.Ot, and t = 0.25) and 
8 (a) and 7 (b) electrons. Lower panels are the NRG results 
with q = 0, iV c = 13 (c) and N c = 14 (d). (e d = -0.5, 
U = 1.0, and V = 0.2.) 



where t is the hopping matrix element in the chain and 
v= 1,2,. ...,7V. 

Using a Lanczos algorithm we first calculate the ground 
state energy and wavevector. In the present case it 

to 



is convenient to use a second Lanczos algorithr 
evaluate directly the impurity spectral function, rather 
than evaluating excited states, their corresponding en- 
ergies and matrix elements M\\>. For a better com- 
parison of the different cases, we always measure the 
frequency from the Fermi energy Ep defined as Ep — 
[E (N e + 1) - E (N e - l)]/2 with E (N e ) the ground 
state energy of the system with N e particles. The spec- 
tral density around the Fermi energy shows a series of 
peaks separated by the characteristic energy A. Each 
peak may be composed by a one or a few delta lines as 
shown in Figs. ||(a) and ||(b). In what follows we dis- 
cuss the low energy structure (uj < A) of the impurity 
spectral function for the case of even and odd number 
of particles corresponding to the at-resonance and off- 
resonance situations respectively. Note that the results 
of Figs. |B|(a) and |5|(b) are not for a half filed system 
and as a consequence the electron-hole symmetry is not 
exactly preserved. 

Even number of particles: The ground state of the sys- 
tem is a singlet and the impurity spectral density shows 
a central peak at oj ~ and two satellites at u> ~ ±A. 
The central peak consists of two lines, one for lo > and 
one for w < corresponding to adding and subtracting 
a particle respectively. The splitting between these two 
lines, that decreases as Vq decreases, is not given directly 
by the confinement energy A but by the energy gained by 
forming the singlet, i.e. is given by the Kondo tempera- 
ture of the finite system. The peaks at u> ~ ±A are also 
made of a couple of delta lines each. The justification of 
this description of a central line and satellites at ±A is 
given by the fact that for smaller hybridizations the two 
central lines approach lu ~ and the impurity spectral 
density reproduces the structure of the unperturbed local 
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density of states (the at-resonance situation of Fig. 0(a)). 
Here we used these parameters for a better comparison; 
results with a smaller hybridization calculated with the 
NRG are shown below. 

Odd number of particles: The ground state of the sys- 
tem is a spin- 1/2 doublet. The spectral densities have 
a low frequency gap with peaks at uj ~ ±A/2 . Each 
of these peaks at u> ~ ±A/2 consist of two delta lines 
that correspond to final states with different total spin: 
when an electron is created or destroyed on the spin-1/2 
ground state, the final state may have total spin zero or 
one. The energy of the final state depends on its total 
spin. Again in this case the general structure is that of 
the underlying local density of states. 

These results should serve us as a guide for the numer- 
ical renormalization group calculation. Since the NRG 
involves some approximations and will be extended to 
more realistic cases, as is the case of a system with a fi- 
nite barrier, it is important to have this exact result as a 
reference calculation. 



B. Numerical Renormalization Group: 
Temperature Results 



Zero 



For a finite system, the renormalization procedure is 
truncated at an iteration N c . In Figs. ||(c) and ||(d) we 
present the NRG results for the impurity spectral den- 
sity of a finite system with an even and an odd num- 
ber of particles respectively. As for the calculation of 
the thermodynamic properties, the system with an even 
(odd) number of particles is evaluated with an odd (even) 
number of shells N c . The spectrum consists of a discrete 
collection of delta lines and, as in the previous section, 
the smoothing is done only for practical purposes in the 
presentation of the data. Our low frequency NRG results 
compare very well with the exact results of a finite sys- 
tem [Figs. |'(a) and |(b)]. A gain, for an even number of 
particles the impurity spectral density consists of a cen- 
tral peak, composed by two delta functions, and the two 
satellites at uj ~ ±A. For an odd number of particles 
the approximate NRG results clearly show the central 
pseudogap and the two satellites at uj ~ ±A/2 with a 
structure that comes from the spin-dependent energy of 
the final state. 

For the remainder, we focus on the more interesting 
case of a system in contact with a macroscopic reservoir. 
This is described with a non-zero parameter a represent- 
ing a finite wall. 

In Fig. ^(a) the impurity spectral density for the at- 
resonance situation is shown for different values of a. For 
a = the isolated cluster results are reproduced with a 
peak at low frequencies. This peak is well separated from 
the Fermi energy due to a large hybridization V used in 
the calculation for practical purposes. In the figure only 
a detail of the low frequency structure is shown and the 
satellite at uj ~ A is not observed. As a approaches 
one, corresponding to an infinite homogeneous system, 
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FIG. 6: Impurity spectral density pi(uj,T = 0) for the at- 
resonance case with N c = 19 (a) and the off-resonance case 
with N c = 20 (b) and different values of a: a — 0.2 (dashed 
line), 0.3, 0.4 , 0.6, 0.75, and 0.9 (thick line). (e d = -0.5, 
U = 1.0, and V = 0.2.) 



the two structures merge into a single Kondo peak cen- 
tered at the Fermi energy. The impurity spectral density 
when the Fermi energy is off-resonant is shown in Fig. 
^(b). For a small a the excitations around uj ~ A/2 
are clearly observed. As a increases the uj ~ A/ 2 struc- 
ture is washed out and simultaneously a narrow Kondo 
resonance develops at the Fermi energy. 

For more realistic parameters (a smaller hybridization 
and consequently a smaller Kondo temperature) the at- 
resonance and off-resonance impurity spectral densities 
are shown in Fig. [?]. These results show that, for the gen- 
eral case of a Kondo impurity in a nanoscopic system with 
an intermediate barrier and the Fermi energy at a reso- 
nant state, we should expect the low frequency spectral 
density pi(uj,T = 0) to present a broad Kondo resonance 
with some structure. In fact, for large barriers (small a), 
the Kondo resonance has a minimum at the Fermi energy. 
For the off-resonant case, we expect the spectral density 
to present the structures at uj ~ ±A/2 separated by a 
pseudogap and - at zero temperature - a central Kondo 
peak. As we show in the next section, a small tempera- 
ture completely destroys this central Kondo peak while 
the broad structures survive up to much higher temper- 
atures. 



C. Phenomenological approach 

Here we briefly discuss the structure of the Kondo res- 
onance in terms of what would be obtained with a simple 
slave boson theory in the saddle point approximation. In 
this theory (U — opj), the low energy impurity Green's 
function is given by 



G(uj) 



ex 



b 2 V 2 g° {uj) 



(9) 



where b 2 is the square of the mean value of the boson 
field, E\ ~ gives the position of the Kondo resonance 
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FIG. 7: Same as in Fig. | with a = 0.45 and V = 0.15. 
a) At-resonance case with N c = 19 and b) Off-resonance case 
with No = 20. 
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FIG. 8: Phenomenological results for the impurity spectral 
densities. (A = 1.0, 7 = 0.1, and b 2 V 2 = 0.02). a) At- 
resonance. b) Off-resonance. 



and g°{oj) is the bare local propagator of conduction elec- 
trons. Here as in the NRG calculation we take the Fermi 
energy equal to zero. In an homogeneous system it is 
a good approximation to take g (u>) — —iirp where p 
is the frequency independent local density of states and 
defining the Kondo temperature T# as the width of the 
Kondo resonance the local propagator can be put as 
G(u) ~ b 2 /(u + iT K ) with b 2 = T K /irV 2 p in agreement 
with the Friedel's sum rule. For our case of a central grain 
weakly coupled with a reservoir, the propagator g°(u)) is 
taken as a sum of poles separated by the characteristic 
energy A and widths 7 

3°H = E ( 10 ) 

*Y UJ - A; + IJ 

with A; = I A or A/ = (/ + 1/2) A for the at-resonance 
and off-resonance cases respectively. The structure of 
the Kondo resonance so obtained is shown in Fig ||. The 
results are in good qualitative agreement with the NRG 
results. In the at-resonance case, depending on the value 
of the parameters, the low energy spectrum consists of 
two peaks with a deep at the Fermi level or a single peak 
with a maximum at the Fermi level (not shown). In the 
off-resonance case, a central peak at the Fermi level is 
obtained. Note however that the area of the central peak 
relative to the area of the satellite peaks obtained is this 
approximation is different from that obtained with the 
NRG. 

D. Finite Temperature Results 

The NRG calculation of the finite temperature spectral 
density p%{u>, T) relies on the same approximations of the 
T = case. The spectral density at a fixed temperature 
T is evaluated as above - using Eq. (§) instead of Eq. (g) - 
if u; > fcflT. To calculate the spectral density at id < fc^T 
a Hamiltonian Hpj with N such that tDjy ~ kgT is used. 

When T£? ~ A and the Fermi energy lies between 
resonances, the off-resonance case, two energy scales are 
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FIG. 9: Impurity spectral density pi(uj,T) for different tem- 
peratures, (a) At-resonance: T — (thick line), T = 0.23A 
(thin line), and T = 0.7 A (dashed line) (b) Off-resonance: 
T = 0.0 (thick line), T = 2.0 x 10 _3 A (thin line), T = 
5.1 x 10~ 2 A (dashed line), and T = 0.5A (dotted line). 

clearly observed in the temperature dependence of the 
susceptibility. In fact as the temperature is lowered they 
determine the onset of the plateau in fc^Tx, that is given 
by A, and the offset at Tjf that is determined by a. 
For the at-resonance case the static magnetic properties 
are dominated by the larger scale A since for fc^T < 
A the magnetic moment is completely screened. The 
spectral densities are also sensitive to the temperature 
and as we will see, there are changes in T) each time 
the temperature approaches each one of the characteristic 
energy scales. 

In Fig. ||(a) the spectral density pi(u>,T) is shown at 
different temperatures for a system with the parameters 
as in Fig. |[ The Kondo peak develops and reaches its 
zero temperature value as soon as the temperature goes 
below A. The parameters are the same as in Fig. || with 
a = 0.35. 

For the off-resonant situation the temperature evolu- 
tion is shown in Fig. ^|(b). In this last case the two 
characteristic temperatures are also reflected in the mag- 
netic susceptibility from which we can extract their val- 
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ues as the onset and offset of the plateau in p? . As can 
be seen from the figure, as the temperature is lowered, 
the broad satellite at cu ~ A/2 develops at a tempera- 
ture corresponding to the onset of the plateau while a 
narrow Kondo structure develops at a temperature cor- 
responding to the offset of the plateau. The first and sec- 
ond temperatures at which the narrow Kondo peak and 
the broad structure disappear indicate the decoupling of 
the impurity spin with the electron spin density at large 
(r > R c ) and short (r < R c ) distances respectively. 

IV. SUMMARY AND DISCUSSION 

We have presented a model for a Kondo impurity in 
a nanoscopic sample or grain coupled with a reservoir 
through a large barrier. The natural energy - or length - 
scale associated to the Kondo effect in a bulk material is 
to be compared with the new scale introduced by the con- 
finement of electrons. These new scales modify the local 
density of states at the impurity site and consequently 
their thermodynamic and spectral properties. 

The local density of states of the host material consists 
of resonant states separated by a characteristic energy 
A and with a width 7. While A is determined by the 
size (R c ) of the grain, the width of the resonances is 
determined by the transparency (a) of the barrier. We 
have considered cases in which the Fermi energy lies at a 
resonance or between two of them. 

The main results of the paper are the frequency and 
temperature dependence of the impurity spectral density 



Pi(uj 7 T). We have performed a detailed analysis of the 
structure of pi(tu,T) for u> around the Fermi energy for 
the case T^? ~ A. This regime with Tj? ~ A is char- 
acterized by a interplay between electron-electron corre- 
lation and confinement effects. The width of the Kondo 
resonance (Tg?) and the characteristic energy (A) that 
defines the structure of the local density of states are of 
the same order and as a consequence the Kondo reso- 
nance shows a superstructure. When the Fermi energy 
lies at a resonance, the low temperature spectral density 
consists of a broad structure around the Fermi level. For 
large barriers, the spectrum presents a two peaked struc- 
ture with a deep at the Fermi level, where the distance 
between peaks is given by the Kondo energy of the finite 
system. As the height of the barrier is decreased, the 
peaks merge into a single broader peak. 

For the off-resonance case Pi(u),T — 0) presents struc- 
tures at w ~ ±A/2 and a central Kondo- like peak at 
lo ~ 0. In temperature dependence two energy scales can 
be distinguished, at which the two types of structures 
disappear. These two energy scales A and T^* are also 
reflected in the magnetic susceptibility. Phenomenologi- 
cal results based on the slave boson mean field theory for 
the impurity spectrum at zero temperature give results in 
qualitative agreement with the ones obtained with NRG. 

Different experiments like the STM conductance for 
impurities in quantum corrals or small terraces or the 
transport through QD in nanoscopic rings could be used 
to test the theory. 

This work was partially supported by the CONICET 
and ANPCYT, Grants N. 02151 and 99 3-6343. 
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